Thermopower of gapped bilayer graphene 
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We calculate thermopower of clean and impure bilayer graphene systems. Opening a band gap 
through the application of an external electric field is shown to greatly enhance the thermopower of 
bilayer graphene, which is more than four times that of the monolayer graphene and gapless bilayer 
graphene at room temperature. The effect of scattering by dilute charged impurities is discussed in 
terms of the self-consistent Born approximation. Temperature dependence of the thermopower is 
also analyzed. 

PACS numbers: 81.05.Ue, 73.50.Lw, 73.50.-h, 72.10.-d 
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I. INTRODUCTION 

Seebeck coefficient, also called as thermopower, mea- 
sures the voltage drop across a material in response to 
a temperature drop. The achievement of large ther- 
mopower is a prerequisite to realistic applications in 
heat to electric energy conversion. Thermopower, among 
other thermoelectric properties, also complements the 
conductivity in elucidating mechanisms dominating the 
transport processespj. Recently, thermopower of mono- 
layer graphene, a peculiar two dimensional electronic sys- 
tem characterized by a Dirac like relativistic dispersion, 
attracts much attention from both experimental [2hJ1 and 
theoretical groups dim. Calculations taking the effect 
of charged impurity scattering into account could explain 
the experimental results very well. The experimentally 
observed deviation from the Mott relationship at low car- 
rier density is interpreted in terms of electron-hole pud- 
dle formationfl'l and also by mixing of valence band and 
conduction band states by impurity scattering 7]. 

Bilayer graphene is another interesting material sys- 
tem displaying many unusual properties. Upon applying 
an external voltage, a semiconducting gap is induced in 
the otherwise zero gap band structure |9l4l4l|. The gap be- 
ing tunable by external potential difference between the 
two layers introduces a new degree of freedom to bilayer 
graphene. Up to now, there is neither experimental nor 
theoretical work on thermopower of bilayer graphene. It 
is the purpose of this work to partially fill this gap by 
theoretically predicting the behavior of thermopower in 
bilayer graphene systems. 

It is established that charged impurity scattering is pri- 
marily responsible for the transport behavior observed 
in monolayer graphene [isl - flTj . For bilayer graphene, 
the prediction in terms of charged impurity scattering 
is shown to be in qualitative agreement with the ex- 
perimental result of the conductivity, and the opening 
of a gap in biased bilayer graphene is proposed to fur- 
ther improve the agreement 18]. In the present work, we 
theoretically study the thermopower of gapped bilayer 
graphene. We treat charged impurity scattering in terms 
of the self-consistent Born approximation (SCBA)0, [l^ 
[2]| . To ensure the applicability of SCBA, we restrict our 
calculations to relatively clean systems with low impu- 



rity concentrations, where the localization effect is not 
severe [22| - [23 |. Thermopower as a function of carrier con- 
centration is mainly calculated at room temperature. We 
also study the temperature dependence of thermopower. 



II. MODEL AND METHOD 

We consider a bilayer graphene system composed 
of two grap hene single layers arranged in the Bernal 
stackingpS'. We start from a tight binding model in- 
corporating nearest neighboring intralayer and intcrlayer 
hopping terms. An on-site potential energy difference be- 
tween the two layers is included to model the effect of an 
external voltage. In the presence of impurity, the Hamil- 
tonian consists of two parts: H = Hq + Himp- Without 
of magnetic field or magnetic impurities, the two spin 
fiavors are degenerate. We ignore the spin degree of free- 
dom here and multiply the results by two for spin depen- 
dent quantities. The free part of the Hamiltonian is then 
written as 



(1) 



in which the vector of fermion creation operators is de- 
fined as = ("ik'^ik'4k'"2k)- "ik ^^d create a 
layer states with wave vector k on the A and B sublat- 
tice, respectively. Up to nearest neighbor hopping, 77o(k) 
is written as[9l. [25l - [28j 
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~t ^ describes the intralayer nearest 
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neighbor hopping with strength t. The three near- 
est neighbor vectors are defined as = (i, ^)a, 

82 = (i,-^)a and <53 ^ (-l,0)a[2|], a=1.42 A is the 
shortest carbon-carbon bond length. t± is the nearest- 
neighbor interlayer hopping energy. In this work, we take 



2 



t=3 eV and t±=0.3 eV. V is the potential energy dif- 
ference between the first and second layers induced by 
a bias voltage. Since for every attainable carrier den- 
sity, it is possible to find a bias voltage to make the 
potential difference between the two layers as V (when 
the gap induced by V is experimentally reachable), we 
would not consider the Coulomb interaction between im- 
balanced electron densities of the two layers and also ne- 
glect the dependence of V on the carrier density n in this 
work[i[l3|2iil|. 

For charged impurities, the impurity scattering part of 
the Hamiltonian is written as0, 0, [lo, l32j 



Vo — 
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Vq is the volume of the system. The charge density op- 



represents the inverse temperature 1/fcsT, with fc^ the 
Boltzmann constant. The value of x thus qtf de- 
pends on both the temperature and the chemical poten- 
tial. 

In order to calculate the thermopower, we should first 
obtain the particle current and heat current operators 
jjv and jg. They are obtained in terms of the continuity 
equation, which for the particle current reads [34, 35] 



p(r)+V-jAr(r)-0. 



(9) 



The momentum space version of the continuity function 
is 



*p(q) = [pi<i):H] = q • jw(q). 



(10) 



A similar relationship holds for the energy density /i_E(q) 
and the energy current j_E(q). With the particle den- 



erator is defined as p(q) = EV-^^k+q- The electron- ^.^^ ^^^^^^^^ ^^^^^^^ ^^^^ = E E («L«",k^ 



impurity scattering amplitude 14 (q) could be written as 
Vi{q)pi{—q), where jOi(— q) and Vi{q) are the Fourier 
components of the impurity density and the electron- 
impurity potential, respectively. For charged irnpurity, 
Vi (q) is taken as of the Thomas- Fermi type[l|, 0, H^l 



k a=l 



^ak^a,k+q), the particle current operator is obtained at 
the q=0 limit as 



jAr(q = 0) =^Vlj?V'k, 



(11) 
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where the matrix is defined as 



e is the effective dielectric constant from lattice and sub- 
strate, e=3 is adopted in this work[3, [13, H^l- d is the 
distance between the impurities and the graphene plane 
and would be set as zero in the present work[l|, Q. qxF 
is the Thomas-Fermi wave number and is obtained from 
the long- wavelength-limit static polarizability of the cor- 
responding noninteracting electron svstempi [20j 



as 



with the static polarizability 
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A factor of '2' comes from the two fold degeneracy in 
spin. The subindex 'c' means retaining only connected 
Feynman diagrams in evaluating the expectation value. 
The particle number operator is defined as 

= E V'l W^k(r) = E ^Ur)^k(r), (7) 

k 
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where (fl. = (c{^, cj^, c^^), with c^^ representing the 
creation operator of the a-th (a=l, 2, 3, 4) eigenstate 
of Ho{k) with eigenenergy denoted as eka- Thus x is 
obtained in terms of the free particle eigenstates as 
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where fkQ=ekQ — P and nF{x)=l/{e^'^ + 1) is the Fermi 
distribution function, fj, is the chemical potential. /3 
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where the velocity is defined as Vk=Vk0(k). The 
momentum space energy density operator could 

be written as[3i| hsiq) = E V'k'^o(k, q)'0k+q + 

fc 

^ Vi{c[')p{c[' + q), where the free part is 
q' 
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(13) 

where 0(k, q) — ((/)(k) -|- (/)(k -I- q))/2. Calculating the 
commutator between /iB(q) and the Hamiltonian, the 
energy current operator could be obtained. The heat 
current operator, defined by the q=0 limit of jqiq) ~ 
j£;(q) - Mjw(q); is written as 



jQ(q = 0) = E^kj2'/'k + -j^ E ^^(q')4j?'''>k+q', 

(14) 



k,q' 



where the free part is written as 
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with d(k) = i[vk(/'*(k) +v*0(k)]. Substituting Vk by 

^(vk+q' + Vk), we obtain j^ '^ from j^. 

We calculate the thermopower in terms of the Kubo's 
formula [H, [s^l with impurity scattering treated to the 
order of self-consistent Born approximation (tI [l9l |32|. [ssj . 
The thermopower is given by 36, 3^, H^] 



Li2 



eTLu 



(16) 



where e is the absolute value of the electron charge. The 
linear response coefhcients Lij are obtained from the cor- 
relation function Cij{iuj) by 



Lij = lim Re£.y + iO"*"). 



(17) 



In the Matsubara notation, the correlation function 
reads Hal 



iT 



{iuJn)dVo Jq 



dTe*'^"^(rj,(T).j,(0)), (18) 



where d = 2 is the dimensionality, /3 — l/kgT, and in- 
dicates an ordering of the current operators with respect 
to the complex time r. cj„ is the bosonic Matsubara 
frequency related with the current operator. 

For charged irnpurities, only intravalley scattering is 
important 0, [l^, [2^ |4l|, so we neglect the intervalley 
scattering processes in this work. Remembering that only 
states close to the chemical potential contribute to the dc 
transport, we could focus our attention at a single valley 
and retain only those low energy states. Here, we would 
focus on the valley around K = (v^, 1) . Similar to 
the monolayer graphene case, this could be achieved by 
introducing an energy cutoff Ec, such that only when 
the smaller positive eigenenergy is less than Ec would 
the states labeled by the corres pon ding wave vector k be 
retained in our calculations [tI, Ti9| . For experimentally 
relevant carrier densities, Eq is no larger than 1 eV. In 
this energy range, 0(k') could be expanded as a polyno- 
mial series in terms of the relative wave vector k = k'— K. 
Here, we retain the expansion to the second order of k as 

3 

(/)(k) = (/)(k + K) = -<^e'('^+^)"'^" 

~ —es^lkyU — ikxtt + -{kyG + ikxtt)'^] (19) 

Test calculations show that the exact dispersion around 
K could be excellently approximated by the above ap- 
proximation up to Ec ~ 1.2 eV, so it is accurate enough 
for our problem. 

In terms of the approximate <^(k), the velocity vector 
satisfies q • (vk -I- Vk+q)/2 = (j){k -|- q) — 0(k). The im- 
purity scattering part of the heat current operator could 
thus be written as ^ ^j(q)jAf(q)- 



Since the particle current operator is the same as that 
without of impurity scattering, the linear response coef- 
ficient Lii could easily be shown to be[3. 
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]Re{Pii(e- iO+,e + iO+) 
(20) 



The kernel is defined as 



Pn {z, z')^—Y, Tr{Gk(^)ri(k, z')G^{z') ■ 

° k 

(21) 

with ri(k) as the vertex function corresponding to the 
wave vector k. 



Taking into account of the relationship '^ '^ — ^j^" 

^Ji^'^i and following the same route as for the single 
orbital model[365, it could be proved that L12 could be 
written into the following form j^, 0, 
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(22) 



The kernel P12 is simply Pn multiplied by the energy 

Pi2(e T «0+, e + ^0+) = ePn{e T «0+, e + ^0+). (23) 

In our calculations, the positive infinitesimal 0^ would 
be replaced by a small positive quantity rj. In this work, 
we adopt 77=1 meV. A smaller rj is verified not to change 
the results to be presented in what follows. 

Under the SCBA, the 4x4 Green's function matrix is 
defined as Gk{z) = (Gliz)-^ - T,i,{z)). The self energy 
is determined by the following self-consistency relation 



^k(^) = ^ E - k')nG°,(z) - Ek'(^)]-\ (24) 



where Gk(z) — [2 + /i — Ho(k)]^^ is the free Green's 
function, fi is the chemical potential determined from 
the free Hamiltonian for a certain carrier density. In this 
work, we would neglect the shift in by the impurity 
potential. This is known as not giving rise to qualita- 
tive changes when the impurity concentration is not very 
highQ- After obtaining the Green's functions, the vertex 
functions are calculated by the following self-consistency 
relation 



ri(k,z,z') 



^|z;,(k-k')pGk'(z)ri(k',z,z')Gk'(^')- 

(25) 



In obtaining the above self-consistent relations, averages 
over impurity configurations have been done under SCBA 
as(37t 



q,q'' 



(26) 



where Ni = riiVo is number of impurities in the system 
under consideration. For a set of k vectors, we can get 
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the Green's functions and the vertex functions, then they 
are used to calculate the kernels Pn and P12, and hence 
the linear response coefficients Ln and L12 could be ob- 
tained. If the full Brillouin zone (BZ) is utilized, the 
number of wave vectors would be too large for a practi- 
cal calculation. However, since only intravalley scattering 
is relevant for charged impurities (tI [T5l-[l7| . we could fo- 
cus on the low energy states within a cutoff energy Ec 
around a single valley. 

III. RESULT AND DISCUSSION 

To see the effect of impurity scattering more clearly, we 
first present the results for clean systems, where the bare 
Green's functions and vertex functions are used. The 
room temperature (300 K) thermopower of clean mono- 
layer and gapless bilayer graphene are presented in Fig 
1(a). Here and later, the abscissa index x represents 
the electron doping averaged to per site, which could 
be controlled by an external gate voltage [2|-|4|. a;=0.001 
amounts to an electron density of 3.82x10^^ cm~^ per 
layer. For the monolayer system [i^, a tight binding 
model up to nearest-neighbor hopping t=3 eV is used. 
Peak values of the two curves are almost the same, with 
that of the monolayer slightly larger. When a more 
realistic parameter t=2.7 eV is used for the monolayer 
graphene [111, the peak position shifts slightly away from 
x—0, while the peak value keeps almost unchanged as 
about ±83 /iV/K. Besides the low carrier density peak, 
thermopower of gapless bilayer graphene shows a second 
smaller peak at a higher carrier concentration. The on- 
set of the second peak, for which the peak position corre- 
sponds to a chemical potential of |/i|~316 meV, is associ- 
ated with the crossing of the chemical potential with the 
lower valence band (band top at -t±=-300 meV) or the 
upper conduction band (band bottom at t±=SOO meV). 

We now explore the effect of opening a band gap in 
the clean bilayer graphene. Figure 1(b) shows the ther- 
mopower of bilayer graphene for a series of potential dif- 
ferences V between the two layers. With the increase of 
V, peak value of the thermopower increases quickly. For 

V = 1 eV (corresponding to an energy gap of approxi- 
mately 288 meV, which is experimentally achievable [l3|). 
peak value of the thermopower is about 412 /itV/K, which 
is more than four times that of the value in monolayer 
graphene and zero gap bilayer graphene. The much 
smaller second peak shifts continuously to larger x as 

V increases and becomes irrelevant as V increases to 1 
eV. Hence, we would concentrate on the region of low 
carrier density in the following. The smallest energy gap 
between the conduction band and the valence band in- 
creases with V as A = [22l . W3 , which is shown 

as an inset of Fig. 1(b). The large magnitude of ther- 
mopower and the tunability of gap make the biased bi- 
layer graphene a more promising candidate for future 
thermoelectric applications as compared with monolayer 
graphene. 
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FIG. 1: (a) Room temperature thermopower as a function 
of carrier density for clean monolayer graphene (MLG) and 
gapless bilayer graphene (BLG). (b) Room temperature ther- 
mopower of BLG for a series of potential energy differences 
V between the two layers. Inset of (b) shows the evolution of 
the global energy gap A as a function of V. 



In the following, we focus on the biased bilayer 
graphene system with V—1 eV. In Fig. 2(a), ther- 
mopower of this system is shown as a function of chemical 
potential (/i) for three different temperatures. When /i is 
at the band edge (^ ±0.144 eV), thermopower of all three 
temperatures are nearly identical to each other. When 
fi is inside of the band gap, |S'| increases as temperature 
decreases. While when ^ lies in the bands, 15*1 increases 
as temperature increases. This is similar to the corre- 
sponding behavior in semiconducting armchair graphene 
nanoribbons[42i |4^. The corresponding dependence of 
X on /i is illustrated in Fig. 2(b). When x is very close 
to zero, decreases as temperature increases. As \x\ 
increases beyond a certain critical value, increases as 
temperature increases. The above complex temperature 
dependence is a direct result of the existence of the band 
edge Van Hove singularities. From Fig. 2(a) and Fig. 
2(b), it is clear that for 100 K, peak value of 5* is achieved 
at a doping very close to zero. As temperature increases, 
the carrier density for the peak increases continuously to 
larger values. 

Now, we begin to study the effect of impurity. In this 
work, charged impurity is considered as the only source 
of scattering for impure bilayer graphene. Thermopower 
of bilayer graphene for a series of different impurity con- 



5 




FIG. 2: (a) Thermopower as a function of chemical potential 
for a biased bilayer graphene, at three different temperatures, 
(b) The corresponding variation of doping concentration x as 
a function of the chemical potential. Inset of (a) shows an 
enlargement of the curves close to the valence band top (~- 
0.144 eV). Lower inset of (b) shows an enlargement of the 
small fj, and small x region. Upper inset of (b) shows an 
enlargement of the region where the temperature dependence 
of for fixed x changes. 



centrations are shown in Fig. 3(a), at room temperature 
for = 1 eV. Curve for clean system is also displayed as 
a reference. In order to ensure that the SCBA is valid, we 
consider only cases with small impurity concentrations. 
Up to a concentration of ni=5x 10^° cm~^ per layer, peak 
value of the thermopower increases continuously with the 
impurity concentration. The peak position remains un- 
changed. In these relatively clean systems, the influence 
of localization [22| - [23 | can be safely neglected. 

Here and later when impurity scattering is considered, 
an energy cutoff of Ec— 0.5 eV is used. Test calcula- 
tions by increasing Ec show no perceivable change in 
the results for both clean and impure systems in the con- 
sidered low carrier density region. For most energies, 
the self-consistency for the Green's functions and vertex 
functions converge within 100 iterations with an accu- 
racy of 10~^ and 10"'' for the modulus of every element, 
respectively. For a symmetric band structure, as is the 
case for both clean monolayer and bilayer graphene, ther- 
mopower is an odd function of the carrier density [11, 0]. 
In the presence of charged impurities, the electron-hole 
symmetry of the band is preserved which could be seen 



FIG. 3: (a) Room temperature thermopower of bilayer 
graphene for a series of different impurity concentrations with 
V—1 eV. (b) The Thomas-Fermi screening wave vector at 300 
K and 50 K for bilayer graphene with V=l eV. The curve for 
50 K is divided by a factor of 2. Insets of (a) and (b) show an 
enlargement of the corresponding low carrier density parts. 



by explicitly calculating the density of states. So we ex- 
pect the relationship S{—x) = —S{x) survives. The full 
curve for ni=5xl0^" cm~^ and T — 300 K is calculated 
explicitly, verifying the above statement. To save com- 
puting time, for all other parameter sets in the presence 
of impurity, thermopower is calculated explicitly for the 
hole doped cases with x < 0. The results for x > are 
obtained through S{x) — ~S{~x). 

The Thomas-Fermi screening wave vector is shown in 
Fig. 3(b) as a function of carrier density at 300 K and 
50 K for V=l eV. At 300 K, as a result of thermal exci- 
tations, we get a finite (/TFa— 0.017 for zero doping. As 
temperature goes down, the zero doping Thomas-Fermi 
screening wave vector decreases gradually. As could be 
seen in the inset of Fig. 3(b), grrix — 0) is very close 
to zero at 50 K. This result is similar to the monolayer 
graphene system and is different from the gapless bilayer 
graphene, for which the zero doping qxF is finite even at 
zero temperature for the nonvanishing density of states 
there [45|. The two peaks in Fig. 3(b) arise from the Van 
Hove singularities near the conduction band bottom and 
the valence band top, and are shifted to larger carrier 
densities at higher temperature by thermal excitations. 

Previous works [ll-H. Fzl. [i^ on monolayer graphene show 
that impurity scattering is essential to reproduce the tem- 
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FIG. 4: Thermopower of gapped bilayer graphene at 300 K 
and 50 K as a function of carrier density, for V=l eV and 
71^=5x10^'' cm-^ 
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perature dependence of the thermopower observed ex- 
perimentaUy. Thermopower of gapped bilayer graphene 
with the same impurity concentration 71^=5x10^'' cm^^ 
are shown as a function of carrier density for room tem- 
perature and 50 K in Fig. 4. Results for x >10~^ are 
readily obtained for both temperatures. However, for 
X at and very close to zero, only results for r=3GG K 
are obtained within our calculation time. This is under- 
stood from Fig. 3(b) as a result of reduced screening at 
a temperature as low as 50 K^t']. As temperature de- 
creases, thermopower is suppressed and peak position of 
thermopower shifts slightly towards zero doping. These 
results are qualitatively very similar to the experimental 
results for monolayer graphene [2 4]. 

A peculiar feature of our results is that the maximum 
value of S at room temperature increases with impurity 
concentration n^, as shown in Fig. 3(a). By increasing 
the energy cutoff Ec and the number of wave vectors in 
the Brillouin zone, we have verified that the above re- 
sult is robust. In order to have a better understanding, 
we show the variation of Ln and L12 for a::=-3xl0^^ in 
Fig. 5(a). It is clear that, both Lu and L12 decrease 
sharply as rii increases. However, the reduction of Lu is 
somewhat larger than that of L12. The increase of room 
temperature S thus comes as a result of the stronger 
dependence of Lu on as compared with L12. The rea- 
son why Lu decreases faster than L12 as rii increases is 
encoded in the integration kernels of the two linear re- 
sponse coefficients, Eq. (21) and Eq. (23). Eq. (23) 
shows that states above and below of the chemical po- 
tential contribute to L12 in opposite sign but of the same 
sign for Lu. Hence Lu is more sensitive to the variation 
of density of states (DOS) in the conduction band and 



FIG. 5: (a) Variation of Lu and L12 as a function of impurity 
concentration for a;=-3xl0~^, normalized by the ni=0 value, 
respectively, (b) Low energy density of states (DOS) for three 
different impurity concentrations. 



valence band due to the presence of impurities which is 
shown in Fig. 5(b) for two impurity concentrations at 
V=l eV and r=300 K. In a previous theoretical work 
on graphene nanoribbonp^ . a similar increase of S with 
defect density is also observed. 

We now try to understand the doping dependence of 
the shape of thermopower in the small carrier density re- 
gion. From the definitions in Eqs. (16), (20) and (22), 
thermopower could be understood as the average value 
of (E — fj.) weighted by the combination of electron group 
velocity (encoded in Vk in the current operator and the 
renormalized current vertex) and DOS |42l . |46| . States 
above and below of the chemical potential contribute in 
the opposite sign to S. At the same time, the factor 
dnp / duj=—nF{i—nF) / is substantially nonzero only 
in an energy range of several ksT centered symmetrically 
around the chemical potential. It is thus easy to under- 
stand that S tends to be larger once the DOS and (or) 
group velocity between states above and below the chem- 
ical potential have a big contrast [l^]. According to this 
picture, when fj. is deep inside of the band (but still in the 
low carrier density region), where the DOS is nearly flat 
according to Fig. 5(b), the thermopower would be small. 
As the chemical potential moves close to the band edge, 
difference between states above and below of the chem- 
ical potential increases, which results in the increase of 
S. At flnite temperature, when the doping is at or very 
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close to the band edge Van Hove singularities, would 
be inside of the gap. In these cases, all states are high 
energy states measured from jj,. Since S would be larger 
once the higher \E — states contribute more to the in- 
tegration, S is expected to continue increase for very low 
carrier densities. As shown in previous works [42i [iq. for 
/i inside of the semiconducting gap, the most significant 
part of S comes from a term of the form 

^-(y-H)/eT, (27) 

which clearly shows the increase of S as /i— >0 (or equiv- 
alently, x— >-0). 

The above picture explains the initial increase of ther- 
mopower as |a;| decreases. However, when the carrier den- 
sity is very close to zero, the chemical potential also lies 
close to zero energy. Hence, both the valence band and 
the conduction band states are present by thermal exci- 
tations. Since the contributions to thermopower from va- 
lence and conduction band states are opposite in sign, the 
thermopower is expected to decrease at a critical carrier 
density characterized by the temperaturefl, H, H^]- The 
critical |a:| below which the magnitude of thermopower 
starts to decrease is thus expected to increase with tem- 
perature. This explains the shift of peak position as ob- 
served in Fig. 2 and Fig. 4. In Ref. 7 on thermopower 
of monolayer graphene, the deviation at low carrier den- 
sity from the higher density Mott's behavior [l], is 
ascribed to impurity scattering mediated coherence be- 
tween the conduction and valence bands. According to 
this mechanism, as impurity concentration Ui increases, 
the above coherence effect should enhance. So a shift 
of peak position with rii is expected. At 50 K and for 
71^=5x10"'^° cm~^, the peak position is shown to shift 
from xc— to xc— 2xl0~^. However, for 300 K, the xc 
show no perceivable variation with Ui up to 71^=5x10^" 
cm^^. A calculation beyond the SCBA is needed to know 
whether or not the peak position for 300 K would shift 
for much larger n^. 

It is also interesting to ask why introducing a gap sig- 
nificantly enhances the thermopower of bilayer graphene. 
Formerly, a 'pudding mold' mechanism [47i] is introduced 
to account for the large thermopower observed in Cobal- 
tates. In that model, a band with a somewhat flat por- 
tion connected to a highly dispersive portion is proposed 
to give high thermopower when the chemical potential 
lies close to the bending point. Band structure of biased 
bilayer graphene is exactly of 'pudding mold' likefgl. [47|. 
So the increase of S with V for a carrier density typically 
of 2;=±0.001, for which the chemical potential lies inside 
of the band, is understood as resulting from the onset 
of the 'pudding mold' mechanism. However, peak value 
of S occurs when situates inside of the band gap. In 
this case, as mentioned above, Smax could be estimated 
by Eq. (27)111 El- For V=l eV, A~288 meV. For 
x=-3x 10"^ and T=300 K, the chemical potential ^~-40 
meV, Eq. (27) gives a value of approximately 347 /iV/K 
which is about 80 percent of the values in Fig. 2(a). So in 



the present case, Smax is set by the energy gap A, which 
increases with V and is bounded by a limit, i^jH, [i^ . 
Hence, the large maximum thermopower in biased bilayer 
graphene is mainly a result of the energy gap. 

Taking into account of the band asymmetry arising 
from the on-site energy difference between the two kinds 
of carbon sublattices, the band gap becomes asymmet- 
ric and the conduction (or the valence, depending on the 
sign of V) band would be more flat[4g|. However, since 
the on-site energy difference derived from experiment is 
only about 0.018 eV|4j], we expect the above effect is 
extremely small and would not change our present re- 
sult much. The interlayer hopping is estimated to be 
in the range of i^^O.3 — 0.4 eV [4^. Our test calcula- 
tion in clean system with t±=QA eV for V—1 eV (with 
other parameters unchanged) gives a peak thermopower 
value of approximately 544 /zV/K at 300 K, which is 132 
/iV/K larger than the result for i^=0.3 eV. So, the ther- 
mopower of gapped bilayer graphene is large regardless 
of the choice of model parameters. 

Applying an electric field to trilayer or other multilayer 
graphene system, a gap could also be induced JBO, 31] . For 
certain parameters, the band structure is also of 'pud- 
ding mold' likejUlsll. It is thus interesting to ask how 
the peak value of thermopower evolves as a function of 
the layer number. In this work, we have only considered 
the effect of charged impurity scatterings. Recently, it 
has been proposed that scattering by short-range disor- 
der may also play an important role in the transport of 
gapless bilayer graphene because the screening of charged 
impurities in zero gap bilayer graphene is much stronger 
than that in monolayer graphene [3) El]- Though we 
believe that the degree of screening in gapped bilayer 
graphene should be much smaller than that in the gapless 
system, it is an interesting question whether the inclu- 
sion of short-range scatterers would change our present 
results much. On the other hand, the regime of large im- 
purity concentration or strong impurity strength deserves 
an explicit study in terms of a less severe approximation 
as compared with SCBA used above. We defer the above 
questions to later studies. 



IV. SUMMARY 

We have theoretically studied the thermopower of bi- 
layer graphene. If a band gap of approximately 288 meV 
is induced in the system by an external bias, the room 
temperature thermopower is greatly enhanced by a factor 
of larger than 4 as compared with that of the monolayer 
graphene and the gapless bilayer graphene. In the pres- 
ence of dilute charged impurities, peak value of the room 
temperature thermopower is shown to increase slightly. 
This behavior is analyzed in terms of the different depen- 
dence of Lii and L12 on the modification of density of 
states by impurity. As temperature decreases, peak po- 
sition of thermopower shifts slightly towards zero carrier 
density in the presence of dilute charged impurities. 
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